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Hybridization effects and multipole orders in Pr skutterudites 

Yoshio Kuramoto, Junya Otsuki, Annamaria Kiss and Hiroaki Kusunose 
Department of Physics, Tohoku University, Sendai 980-8578, Japan 

Theoretical account is given of 4f-electron dynamics and multipole orders in Pr skut- 
terudites with particular attention to (i) mechanism of the crystalline electric field (CEF) 
splitting leading to a pseudo-quartet ground state; (ii) Kondo effect due to exchange interac- 
tions involving the pseudo-quartet; (iii) multipole orders in the lattice of the pseudo-quartet 
in magnetic field. Competition between the point-charge interaction and hybridization be- 
tween 4f and conduction electrons is identified as the key for controlling the CEF splitting. 
It is found that one of two pseudo-spins forming the pseudo-quartet has a ferromagnetic ex- 
change, while the other has an antiferromagnetic exchange with conduction electrons. The 
Kondo effect is clearly seen in the resistivity calculated by the NCA, provided the low-lying 
triplet above the singlet is mainly composed of the /Vtype wave functions. If the weight of 
the /5-type is large in the triplet, the Kondo effect does not appear. This difference caused 
by the nature of the triplet explains the presence of the Kondo effect in PrFe4Pi2, and 
its absence in PrOs4Sbi2. By taking the minimal model with antiferro-quadrupole (AFQ) 
and ferro-type intersite interactions for dipoles and octupoles between nearest-neighbors, the 
mean-field theory reproduces the overall feature of the multiple ordered phases in PrFe4Pi2- 
The AFQ order with the /3-type symmetry is found to be stable only as a mixture of 0% 
and O2 components. 

§1. Introduction 

Novel many-body phenomena in solids often appear with the support of a specific 
structural feature of the system. One of the typical examples is the high-temperature 
superconductivity in cuprates which is supported by the nearly two-dimensional 
square lattice. In this paper we discuss another interesting example of specific struc- 
tural systems called filled skutterudites. A large number of compounds in this cat- 
egory are characterized by a cage-like structure of ligands in which rare-earth ions 
are loosely bound. In filled Pr skutterudite compounds, many anomalous proper- 
ties have recently been discovered such as the heavy-fermion superconductivity and 
high-field ordered phase in PrOs4Sbi2, 1 ^' 2 ^ antiferro-quadrupole order in PrFe4Pi2, 3 ^ 
and metal-insulator and structural phase transitions in PrRu4Pi2- 4 ^' 5 ^ While neu- 
tron scattering has observed clear CEF transitions in PrOs4Sbi2, 6 ^ 9 ^ only broad 
quasi-elastic features are visible in PrFe4Pi2 above the temperature of quadrupole 
order. 10 ) On the other hand, intriguing temperature dependence of CEF levels has 
been observed in PrRu4Pi2 below the metal-insulator transition. n ) Recently, a new 
phase has been found in PrFe4Pi2, 12 ^ which appears only in high magnetic field along 
the (111) direction. For proper understanding of these phenomena, both on-site and 
intersite interaction effects of 4/ electrons should be taken into account. The on-site 
hybridization with conduction electrons gives rise to CEF splittings in the second 
order and to the Kondo effect in higher orders. The intersite interactions give rise 
to various electronic orders. 

The purpose of the present paper is to give a comprehensive, though still crude, 
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account of rich dynamical and ordering phenomena in Pr skutterudites. Particu- 
lar attention is paid to the following aspects: (i) mechanism of crystalline electric 
field (CEF) splitting leading to a low-lying pseudo-quartet; (ii) Kondo effect due 
to exchange interactions involving the pseudo-quartet; (hi) multiple orders of the 
pseudo-quartet lattice in magnetic field. Most of the theoretical results presented 
in this paper have been originally reported in refs.l3)-15). Here we organize these 
results toward a comprehensive picture. 

§2. Mechanism of CEF splittings in Pr skutterudites 

2.1. CEF parameters under the tetrahedral symmetry 

Filled skutterudites RT4X12 form a bcc structure with the space group Im3. In 
each unit cell, rare-earth ion R is surrounded by 8 transition metal ions T which form 
a cube, and 12 pnictogens X which form an icosahedron deformed slightly from the 
regular one. Figure ^ shows the cage structure with R in the center. The R site has 




Fig. 1. The cage structure forming a unit cell of filled skutterudites. . 

the local symmetry Th which has no four- fold rotation axis. 16 ) In this symmetry, the 
4/ 2 Hund's-rule ground states 3 #4 of Pr 3+ split into a singlet r±, a non-magnetic 

doublet I23, and two I4 triplets. Of these, two I4 triplets are written as and 

(2) 

r} , which are linear combinations of triplets i~4 and I5 in Oh- Under this crystal 
symmetry, the CEF potential is written as 



^cef = A 4 [Ol + 50|] + A c 6 [O u 6 - 210|] + A\[Ol - 01} 
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in the standard notation. 16 ) The term yO\ mixes the I4 and i~5 triplet states in the 
point group Oh- 

There are two main sources for the CEF splitting: the Coulomb potential from 
ligands, and hybridization between 4/ electrons and ligands. Let us begin with the 
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Coulomb potential. In the point charge model, CEF coefficients ^4,^4g and Aq in 
eq. (j2-ip are determined by coordination of the charge and the radial extension of the 
4/ wave function. By interpolating the data of ref.18) for the radial extension, we 
adopt the values (r 4 ) = 3.4ag and (r 6 ) = 19a B with ob the Bohr radius. Explicit 
results have been obtained in ref.13) taking the effective charge Z% of a transition 
ion as the parameter. Here we require the charge neutrality with trivalent Pr as 

3 + 4Z t + 12Z P = 0, 

where Z p is the effective charge of a pnictogen ion. The singlet i"i is stabilized 
for Z% > 0. The stability of the singlet for positive Z t implies that the maximum 
energy gain is provided by a more isotropic electron distribution associated with 7\. 
The CEF potential becomes almost isotropic for Z t ~ —0.25. The quasi-isotropy 
is due to high coordination numbers of pnictogens and transition ions, which tend 
to compensate the anisotropy. For negative charge of transition ions, the electron 
distribution tends to be anisotropic to avoid the Coulomb repulsion. In fact, the 

(2) 

triplet r} is stabilized for Z t < —0.25. To estimate the CEF potential, we use the 
lattice parameter in PrC^Sb]^, 17 -* which gives the Pr-T distance d% = 4.03A with 
T=Sb, the Pr-X distance d p = 3A8A with X=Os and the X-Pr-X vertical angle 
28q = 49.2°. These data will be used in the next subsection to derive the results 
presented in FigEJ 

2.2. CEF splitting by p-f hybridization 

Another important mechanism for CEF splittings is the covalency effect, or hy- 
bridization between localized and ligand orbitals. According to band calculation, 19 ) 
the conduction band striding the Fermi level is formed mostly by the molecular or- 
bital a u formed by 12 pnictogens around each Pr. The other relevant molecular 
orbital t u form two energy bands around a few eV above the Fermi level, and three a 
few eV below. We neglect contributions from t u bands to CEF splittings because of 
the larger excitation energy. Therefore we concentrate on hybridization of the form 

H hyh = V 2u Y,f^a + h.c, (2-2) 

where c a annihilates a conduction electron in the Wannier orbital with the a u symme- 
try at the origin, and f& creates an 4/ electron with the same (a u ) orbital symmetry. 
We take the hybridization parameter Vm real. Hereafter we adopt the Mulliken 
notation such as a u for orbital symmetry, and the Bethe notation such as J\ for a 
double-group representation with spin-orbit coupling. In the second-order perturba- 
tion theory, the effective interaction is given by 

H e $ = PHhyb-^ TrQH^yhP, (2-3) 

h, — Uq 

where P is the projection operator onto 4/ 2 states, and Q = 1 — P. We first deal 
with such part of the second-order hybridization that is diagonal with respect to the 
conduction states. Diagonalization of i/ e ff with this constraint gives the CEF wave 
functions and their energies. 
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We assume the following intermediate states: (i) 4/ 1 and an extra electron in 
vacant states, and (ii) 4/ 3 and extra hole in filled states. For simplicity, we neglect the 
multiplet splittings in both kinds of intermediate states. Energy shifts in respective 
cases are given by diagonalization of the following (2 J + l)-dimensional matrices: 

AE-(M>,M) = -(1 - n 2u ) V -f "£(M'\tif a \M), (2-4) 



A\ + e 2u 

2 

lu 

A 3 + NJ 



AE+(M',M) = -n 2u V j u - ^2(M'\f a fl\M), (2-5) 



where A 2 ±\ (> 0) are excitation energies to 4/ 2±1 , and n 2u is filling of the pnictogen 
state a 2u per spin. By adding contributions from 4/ 1 and 4/ 3 intermediate states, 
we obtain 



AE(M',M) = — V 2 U 2n 2u 5 M ,M' 

^3 + \ e 2u\ 



ri2u 1 - n 2u 



A3 + \e 2u I A x + £2u 



J2(M'\tifa\M), (2-6) 



where the first diagonal term does not contribute to the CEF splitting. The matrix 
elements in eq. (|2-6|) can be derived by explicitly using the 4/ 1 intermediate states as 

{M'\ft\m ims ) = -y/2 Yl (JM'\LM L SM S ) 

x(LM L \lmilmi){SM s \8m B 8a)(lrn! l \ir2), (2-7) 

where m, mi denote azimuthal quantum numbers of I = 3, and m s denotes that of 
s = 1/2. The last factor in ea. (|2-7j) describes the projection to the a2 U -type orbital 
of 4/ electrons. We obtain (lmi\ir 2 ) = ±l/y/2 for mi = ±2, and zero otherwise. 

We note that 4/ 1 and 4/ 3 intermediate states give opposite contributions to 
the level splitting. Hence the sequence of CEF levels is determined by competition 
between both intermediate states. 21 ) In ea. (|2-6jl . the competition appears in the 
factor with n 2u . We have also done calculation where only Hund's-rule ground 
states in 4/ 2=tl are considered as intermediate states. The tendency of opposite 
contributions from 4/ 1 and 4/ 3 states persists, although the relative weight is not 
completely the same. 

The hybridization is parameterized by the Slater-Koster parameters (pfcr) and 
(pf 71 )- It turns out hybridization with a u comes only from (p/tt), which gives the 
interaction parameter 

V 2u = -\/3O(p/7r)sin20o. (2-8) 

With only the a u band taken into account, eq.(J2H)J) shows that the CEF level sequence 
is determined by two parameters; occupation of the a u band and A3/A1. To a 
good approximation we can assume the half-filled a u band, and e 2u = for the 
band. According to available information on PrP, 21 ) we estimate the ratio roughly 
as A3/A1 = 0.6. Hence the level structure due to hybridization alone is more 
influenced by 4/ 3 intermediate states than by 4/ 1 ones. 
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2.3. Combination of Coulomb interaction and hybridization 

We combine both contributions to CEF splittings by point-charge interaction 
and by hybridization. Figure El shows computed results as a function of strength of 
hybridization, (p/7r) 2 /Z\_ where we have introduced 

= l/As-V^i- 

For the point charge parameters we tentatively take Zt = 2 and Z p = —11/12. We 
note that in Mossbauer experiment on PrFe4Pi2, 22 ^ Fe is reported to be trivalent. 
On the other hand, the energy-band calculation suggests a smaller valency for Fe. 19 ) 
The choice Zt = 2 is a compromise between these conflicting results. Overall CEF 
splittings become larger with larger Z t > 0. It should be emphasized that the the 
closely lying singlet and triplet are realized only with simultaneous action of point- 
charge and hybridization interactions. 

The lowest triplet It; which is either or , can be written in terms of Oh 
triplets i~4 and I5 as 

\r u m) = y/w\r±,m) + \f\ - w\r 5 ,m), (2-9) 

where m = ±, specifies a component, and w with < w < 1 gives the weight of i~4 
states. We have chosen the phase of wave functions so that positive coefficients give 
the lower triplet. As the parameter (p/7r) 2 /Z\_ increases, the state is stabilized 
and eventually crosses the i~i level. After the crossing, the ground state acquires 
the local moment. In the case of PrOs4Sbi2, the parameter (pfir) 2 /A^ = 190K is 
consistent with the combination of A\ ~ 5eV, A3 ~ 3eV and {pfir) — 0.35eV. This 
value of (pfir) is about 50% larger than that estimated for PrP. 21 ) Since the distance 
dp ~ 3.5A in PrOs4Sbi2 is smaller than the Pr-P distance 4.2A in PrP, the larger 
magnitude is reasonable. In the case of PrFe4Pi2, the lattice constant is 7.81A as 
compared with 9.3A of PrOs4Sbi2- Hence (pfir) in PrFe4Pi2 should be even larger. 

The level repulsion between and around (pfir) 2 /A- = 110K is due to 
mixing of wave functions, and is a characteristic feature in the point group T^. In 
contrast, a level crossing occurs in the Oh group since there is no mixing between 
i~4 and I5. Possible range of the effective hybridization in some representative Pr 
skutterudites is schematically shown by the lines (a)-(d) in FigEJ In PrOs4Sbi2, the 
triplet forming pseudo-quartet is almost of the i~5-type as shown by (c). On the 
other hand, the CEF levels in the high-temperature phase of PrRu4Pi2 corresponds 
to (a), which bifurcates into (b) and (d) in the low-temperature phase. Namely 
in one sublattice called Prl, the triplet changes the character slightly, while in the 
other sublattice called Pr2, the triplet becomes the ground state. Below the phase 
transition, the system behaves like an insulator. It may be possible to describe 
the transition as the antiferro-hexadecapole order. Explicit results of more detailed 
study will be reported elsewhere. Finally, the case of PrFe4Pi2 should be mentioned. 
As discussed later, the CEF levels seem to correspond to the range beyond the left 
end of Fig|21 At first sight, the larger hybridization in PrFe4Pi2 seems incompatible 
with this leftward shift from PrOs4Sbi2- However, since the smaller lattice constant 
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Fig. 2. CEF level structures derived from hybridization and point charge potential as a function 
of {pfirf/A-. The level sequence qualitatively corresponds to: (a) PrRu4Pi2 in the high- 
temperature phase; (b) Prl site in PrRu4Pi2 in the low-temperature phase; (c) PrOs4Sbi2; (d) 
Pr2 site in PrRu4Pi2 in the low-temperature phase. See text for details. 



in PrFe4Pi2 favors the 4_f 1 excited states rather than the 4/ 3 states, the parameter 
A\ becomes smaller, leading to larger Z\_ and smaller (p/7r) 2 /Z\_. 

§3. Exchange interaction in the pseudo-quartet 

3.1. Second- order exchange interaction 

Now we deal with such components of H c r given by ea. (|2-3|) . that give rise 
to exchange-type interaction between 4/ and conduction electrons. We restrict to 
the singlet and the lowest triplet for the CEF states for 4/ states, and neglect the 
multiplet splittings in the intermediate states. We obtain for 4/ 1 as intermediate 
states, 

1 MM 1 era' 

where M denotes azimuthal quantum number of J = 4, and we have introduced the 
notation 

A(MM'; go-') = {M\flf a \M'), (3-2) 

for the matrix element. In the case of 4/ 3 intermediate states without multiplet 
splittings, the effective Hamiltonian is given by 

#cff[4/ 3 ] = E E^'( MM ';^') i m x m 'i 4<v, (3-3) 

3 MM' cro-' 
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where A'(MM';aa') = (M\f a fl,\M'). Because of completeness of the 4/ 1 and 4/ 3 
intermediate states, we obtain the relation 

A(MM'; go-') = 5mmiS„v> - A' (MM'; aa'), (3-4) 

where the first diagonal term is irrelevant. Thus consideration of 4/ 3 intermediate 
states in addition to 4/ 1 is accomplished by the replacement 

Vl/A, - VKl/A, + I/A,) = Vl/A, 

in eq. <|3-1|) . Note that 4/ 1 and 4/ 3 states contribute additively in contrast with the 
case of CEF splittings. Although the bare diagonal terms are already included as 
CEF splittings, higher order potential scattering terms contribute to renormalization 
of CEF levels as well as damping of states. 

3.2. Symmetry analysis of the pseudo- quartet 

We are interested in the matrix element between the singlet and the triplet in 
the 7}j group. As we have seen the triplet is a linear combination of A and A in 
the Oh group. Hence it is convenient to analyze the symmetry properties of triplets 
in the cubic case. The direct products of pseudo-quartets can be decomposed as 

(A © A) © (A © A) = A + 2A + (A + A + A + A), (3-5) 
(A © A) © (A © A) = A + 2A + (A + A + A + A), (3-6) 

where A has the same symmetry as the magnetic moment. The product of A © A 
produces newly two A representations, of which one is time-reversal odd and the 
other is even. The odd representation represents the magnetic moment, while the 
even one the hexadecapole moment. On the other hand, the pseudo-quartet A © A 
does not produce a new A representation. Physically, this means that A as the first 
excited states gives rise to a van Vleck term in the magnetic susceptibility, while A 
does not. 

It is convenient to introduce the effective moment operators in the pseudo- 
quartet in A. In the case of w = (pure A triplet), the magnetic moment within 
A is the only relevant quantity. The vector operator in this case is written as X t . 
On the other hand, another vector operator X s is necessary to describe the mag- 
netic moment of van Vleck type. The operator X 1 act on the triplet states, and X s 
connect the singlet and the triplet as shown in Fig. |3J 

In terms of the four-dimensional basis set | A), I A> +), I A; 0), | A, — ), Each of 
the operators Xj,Xj with i, j = x,y,z can be represented by a 4 x 4 matrix. Alter- 
natively, one can use the pseudo-spin representation as in ref.23): 

X t = S\ + S2, X b = S\ — S2, 

where S\ and S2 are spin 1/2 operators. The pseudo-spin representation naturally 
leads to commutation rules among X 1 and X s . 

We project the basis set \M) to the pseudo-quartet. In terms of two vector 
operators X and X s , we then obtain the following form for the effective interaction 



8 



Y.Kuramoto, J.Otsuki, A. Kiss and H.Kusunose 



Ft) 



X 1 



X s 



Fig. 3. Effective moment operators in the CEF singlet-triplet system. 



within the pseudo-quartet: 

H cxc = (j t X* + / S X S ) • s c , (3-7) 
where ^ 

The coupling constants are given by I a = V2 U a a /A for a = s,t with 

10 + 88w 4 [5w . . 

at = -^155-' as = -33VT- (3 ' 8) 

In deriving eq. (j3-8j) . we have used eq. ([2-7[l . Both a t and a s are negative, which means 
ferromagnetic exchange between conduction and 4/ moments. The origin of the fer- 
romagnetic sign for at is traced to the Hund rule involving the spin-orbit interaction; 
the dominant orbital moment is pointing oppositely from the spin moment. In other 
words, the spin exchange is antiferromagnetic as is usual for a hybridization induced 
exchange. The sign of a s , on the other hand, is not physically meaningful. 
The exchange interaction in terms of pseudo-spins is given by 

H cxc = (JiSi + J2S2) ■ s c , (3-9) 

where J\ = 1% + -^s and J2 = h — Is- Figure 21 shows the exchange interactions J\ and 
J2 of pseudo-spins in units of V^/A. It should be noticed that J2 becomes positive 
for w > 0.00324, and represents the antiferromagnetic exchange. The emergence 
of antiferromagnetic exchange is due to the particular CEF level structure in Pr 
skutterudites. The antiferromagnetic exchange is almost negligible in the pure i~5 
case (w = 0) , and becomes an order of magnitude larger as w increases toward unity, 
i.e., toward pure I4. It is likely that the Kondo-type behavior seen in PrFe4Pi2 
originates from a large value of w together with a small singlet-triplet splitting. We 
shall next address to this problem. 



§4. Kondo effect in the pseudo-quartet 



4.1. Application of the NCA to 4/ 2 CEF states 

We proceed to derive dynamics of the singlet-triplet system by adopting the 
non-crossing approximation (NCA). 24 )' 25 ) In applying the NCA to systems with ex- 
change interactions, we introduce a fictitious 4/ 1 intermediate state with negligible 
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Fig. 4. Coefficients of pseudo-spins in the effective exchange. 



population. 26 ^' 27 ) We can alternatively introduce a fictitious 4/ 3 intermediate state 
instead of the 4/ 1 state with the same results. The NCA utilizes the resolvent R s (z) 
for the singlet state and Rt(z) for the triplet. The effects of the interactions are 
taken into account in terms of the self-energy S a {z). Each resolvent R a {z) is given 
by 



Ra(z) 



E a (z)Y 



(4-1) 



where e a with a = s,t denotes CEF levels of 4/ 2 states, i.e., Aqef = (-t — without 
hybridization. The NCA determines U a (z) in a self-consistent fashion. 

There also appear renormalized exchanges h{z) and I s (z), which are modified 
from the bare ones It and I s in ea. (|3-7|) . In addition, an effective potential K t {z) 
for the triplet and K s (z) for the singlet are generated by higher-order exchange 
scatterings. The self-energy is given in terms of the effective potential by the NCA 
integral equation: 



U a (z) 



dep c (e)[l- f(e)]K a (z-e), 



(4-2) 



where p c (e) is the density of states of the conduction band, and /(e) is the Fermi 
function with the chemical potential being 0. In order to derive equations for the 
renormalized interactions I a (z) and K a {z), we divide the operator X s into two parts: 



X s = P t X s P s + P s X s P t , 



(4-3) 



where P a is the projection operator to the state a = s,t. The first term operates on 
the singlet state, and second one to the triplet states. Correspondingly, we define 



effective interactions Ig (z) and I^'^) for each part. The simultaneous equations 



r(st) 



10 



Y.Kuramoto, J.Otsuki, A. Kiss and H.Kusunose 



for the effective interactions are given in the matrix form by 



/ h \ 

(ts) 

3 

;(st) 



I, 



h 

Is 



V o / 



/ W2 

/tilt 

i s n t o 
i t n t /2 o 
V o 



I s i7 t 3/4 

where the auxiliary quantity II a (z) is introduced by 



i s n s /2 i t n t 

o o 

o i s ii t 

/ s i/ s /4 





\ / It \ 

i s n s 
o 





(ts) 

(St) 



dep c (e)f(e)R a (z + e). 



Is 

K t 

V / 

(4-4) 



(4-5) 



Solving eq. (|4-4|) for the renormalized interactions, we obtain 

4(2/ t - I s 2 iT s ) 



(2 - /tilt) (4 + 4/ti/t - 3/ s 2 i/ t i/ s ) ' 
^ 



J — f(ts) _ f(st) _ 

s ~ 8 s (4 + 4/ti/t -3/ s 2 i/ t i/ s )' 



K s 



4/ t 2 i/ t + 2/ s 2 iZ s - 3/ t / b 2 i/ti/s 
(2 - / t i/ t )(4 + 4/ti/t - 3/ s 2 i/ti/s) ' 

3/ 2 i/ t 
4 + 4/ t i/ t -3/ 2 i/ t i/ s ' 



(4-6) 



The effective interactions I a (z) are required for the dynamical magnetic susceptibil- 
ity. This quantity is computed explicitly in ref.14) with due account of the so-called 
vertex correction. 29 ) In this paper, we only deal with the electrical resistivity for 
which the vertex correction is not necessary. 

4.2. Temperature dependence of resistivity 

We discuss how the CEF singlet together with the exchange interactions influ- 
ences the resistivity. First we note that the CEF splitting is renormalized signifi- 
cantly by fourth-order effects of hybridization, i.e., the second order in the exchange 
interaction. By perturbation theory with respect to / S) t we obtain the second-order 
shift 



A 



(2) 

CEF 



-D P l{ll 

,(2) 



-Dp 2 c JxJ 2 , 



(4-7) 



which also appears as the difference E^ z> (z) — U^\z) of second-order self-energies. 
The renormalization depends hardly on temperature and energy as long as T,\z\ <C 
D. Since we have J1J2 < 0, the second order exchange stabilizes the CEF singlet. 
On the other hand, higher order exchanges cause the Kondo effect. The competi- 
tion between the Kondo effect and the CEF effect depends on their characteristic 
energy scales. In computing dynamical quantities, we take J2 = to minimize renor- 
malization of the CEF splitting. Even with J2 = 0, the pseudo-spin S2 interacts 
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with conduction electrons indirectly through the other pseudo-spin S%. With the 
additional condition Z\cef = 0, the indirect coupling disappears and the singlet- 
triplet Kondo model is reduced to the Kondo model with the pseudo-spin S2 being 
decoupled. 

In the NCA for the exchange interactions, the impurity T-matrix is given by 



Z f Jc 2vrz ^ 



(z + ie n ), (4-8) 



where Zf is the partition function of 4/ electrons with account of interactions, e n = 
(2n + 1)ttT is the Matsubara frequency of fermions, and the contour C encircles all 
singularities of the integrand counter-clockwise. We utilize the two kinds of spectral 
intensities: 25 ^' 28 ) 

n a (uj) = --lmR a (uj + iS), (4-9) 

7T 

= ZjV^TfcH, (4-10) 

where 5 is positive infinitesimal. In a similar fashion, we further define the spec- 
tral intensities t]2^\uj) and ^\u>) for the effective potential K a (z), and rja\uj), 
£,a\u>) for the renormalized exchange I a (z). Actually the spectral functions £ a (w) 
are computed by another set of equations to avoid divergent Boltzmann factor at low 
temperature. 28 - 1 Performing analytic continuation ie n — > uj + i5 to real frequencies 
in eq. (|4-8|) . we obtain 

-hriI{u + i5)=Y J I + u,) + *#°(e)£«(e + u)}. (4-11) 

a J 

The electrical conductivity cr(T) is derived from ImT(w) by 



oc 



a{T) = AI (4-12) 



OO 



where A is a constant and /(e) is the Fermi distribution function. We take a rect- 
angular model with the band width 2D for conduction electrons 

p c (e) = 9(D-\e\)/(2D), (4-13) 

where 6(e) is a step function. We take D = 10 K and use Kelvin as the unit of 
energy. The cut-off energy D is important in determining the Kondo energy scale, 
but otherwise does not enter the physics as long as it is much larger than other 
parameters. Figure [5] shows temperature dependence of the electrical resistivity 
p(T) = l/a(T) with J\p c = 0.2 and J2 = 0. In the quartet case Z\cef = 0, p(T) 
continues to increase as the temperature T decreases. On the other hand, with the 
CEF singlet ground state the enhancement of p(T) is suppressed at temperatures 
below about the third of Z\cef- The suppression is due to a pseudo-gap of magnitude 
~ 2Z\cef in IniT(w) around the Fermi level. Note that the temperature at the peak 
of p(T) is substantially smaller than ZicEF- Namely the Kondo effect persists even 
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though the CEF triplet is thermally depopulated. This is understandable because 
the van Vleck-type polarization is responsible for the Kondo effect here. Although 
the Kondo temperature is difficult to be defined unambiguously in this system, the 
magnetic relaxation rate at low T is only a few K in the case of ZicEF = 0. 14 ) 





A CEF = 0K 




A CEF = 20K 




A CE F = 40K " 









1 10 100 

Temperature [K] 

Fig. 5. The electrical resistivity as a lunction ol temperature for selected values ol /Acef- 



§5. Multiple orders in PrFe4Pi2 

5.1. Splitting of the CEF triplet by magnetic fields 

We are now going to discuss the phase diagram in PrFe4Pi2 as a representative 
of the multiple ordered phases in magnetic field. Let us compare the situation with 
PrOs4Sbi2 where an AFQ order has been found only under magnetic field. In the 
case of PrOs4Sbi2, the CEF triplet consists mostly of J5, and the AFQ order also has 
the /5-type symmetry. In this case the splitting of the i~5 triplet under the molecular 
field is symmetric about the origin. Therefore, the quadrupole order parameters in 
A and B sublattices have the same magnitude. In contrast, the splitting of the J4 
triplet under the molecular field is asymmetric about the origin in general. This 
feature brings about interesting consequences in the AFQ order in PrFe4Pi2, which 
will be described below. 

We begin with analysis of the CEF energy spectrum in magnetic fields taking 
the following single-site Hamiltonian: 

H ss = Vcef — S/^bH • J, (5-1) 

where the second term represents the Zeeman energy with the angular momentum 
J, and the first term is the CEF potential given by ea. (|2-l|) . The relevant parameter 
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range for PrFe4Pi2 is Wx < 0. 13 ) Then we obtain in the case of x > and W < 0, 

W = 1 -(l + , 3 + 2 * V (5-2) 

2^ 7(3 + 2x) 2 + 10081/ 2 / 

and the CEF levels are given by 

E(r t ) = 2W \x - 4 + 2^(3 + 2x) 2 + 1008y 2 
E(r x ) = W (108x - 80) . (5-3) 

The parameter y can make the T\ and 7^ states closely located to each other. 
For example, the choice x = 0.98 and y = —0.197 gives the vanishing singlet-triplet 
splitting. We recall that a constraint |a;| < 1 is required in the parameterization of 
ea. (|2-l[) . and note that PrFe4Pi2 has the parameter is in the range \x\ ~ 1. Because 
of these reasons, there appears apparent discontinuity in the CEF parameters as the 
CEF potential changes continuously. 13 ) In the case of x < and W > 0, the weight 
w and the CEF levels can be obtained by the following replacement in eas. ()5-2(l and 

x -» x/(2x + 1), y -> y/{2x + 1), W -> {2x + l)W. (5-4) 

It is seen that Wx and Wy remains the same after the replacement. In the special 
case of x = —1, the replacement above is just the reversal of signs for x, y, W. 

Now we derive the eigenvalues of the 9x9 matrix for H ss as a function of 
magnetic field. Figure |H1 shows the result for magnetic field along (111) with the 
singlet-triplet splitting Z\cef = 2K. The most important feature is the level crossing 




near H ~ 30T, which was first noticed by Tayama et al. 12 - 1 The level crossing occurs 
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because of the level repulsion between the states originating from the doublet i~23 
and the triplet i^i . Hence, if one neglects i~23 from the start, the level crossing does 
not occur. The (111) direction of the magnetic field lowers the crystal symmetry to 
trigonal. This symmetry is essential in bringing about the level crossing of the lowest 
two levels. For comparison, we have also derived the CEF spectrum for directions 
(100) and (110) of the magnetic field. Figure [7| shows the results. It is evident that 
the lowest two levels never cross in these directions of the field. We shall discuss 
later that the level crossing for (111) is responsible for the ordered phase which 
appears only for this field direction, and for the anomalous angle dependence of the 
resistivity. 15 ^' 30 ) 




5.2. Instability of the paramagnetic phase 

We now consider the simplest intersite interaction that can reproduce the phase 
diagram of PrFe4Pi2 at least qualitatively. For this purpose we introduce notations 
for the quadrupole operators of the /3-type, namely i~23 in the Th group, and such 
octupoles that reduce to r$ u in the cubic symmetry. These are given by 31 ) 

2 = (0° 2 , Of) = (^(2J Z 2 - Jl ~ J 2 y ), Jl ~ 4) > ( 5 - 5 ) 

(T x >1y ,T Z ) = (^J x Jy — JgJx, JyJz ~ J%Jyi JzJx ~ JyJzj ' (^'^) 



where bars on the products represent normalized symmetrization, e.g., J x Jy = 
{JxJy + JyJxJy + JyJx)/3. We neglect the other quadrupole operators of the /5-type 
in the intersite interaction. This is because the induced antiferromagnetic moment 
in in PrFe4Pi2 is parallel to the magnetic field, 37 ) which suggests that the AFQ order 
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is of i~3-type. Then we take the following model for the intersite interaction: 

Hint = (^dipJi • Jj + 5quad0 2 , i • Q 2 j + 5oct Tf • Tf ) , (5-7) 

where {ij} means the nearest-neighbor pair, and the coupling constants have signs 
5quad > 0, g oct < and g^ip < 0. Namely we take antiferro-type quadrupolar and 
ferro-type dipolar and octupolar interactions. The presence of the ferromagnetic in- 
teraction is suggested by the temperature dependence of the magnetic susceptibility, 
and by the ferromagnetism occurring with La substitution. 32 ) The dipole and I3 
quadrupolar interactions are consistent with the invariant form allowed by symme- 
try of the bcc lattice. 33 ^ For simplicity we neglect the mixing term Jj • T| u , although 

the mixing is allowed by the Th symmetry. 33 ) 

We take the model given by H ss + Hi nt , and will obtain the phase boundary in 
the plane of magnetic field and temperature. The nine crystal field levels labeled 
by 1 1) are derived for different values of magnetic field by solving H ss \l) = Ei\l). In 
deriving the phase boundary, we keep only the lowest six levels for simplicity of the 
numerical calculation. As seen from Fig|f3 the three higher levels neglected here 
repel with lower levels, and go up in energy for H > 20T. Thus their neglect will not 
influence the low-temperature behavior significantly. 

From the six levels labelled by \l) with I = 1,...,6, we have 35 (= 6 x 6 — 1) 
independent pairs k, I describing multipole operators at each site. Among these, we 
have 15 symmetric combinations of the pairs, and arrange them as increasing order 
of Wk + l with k < I. Then we define X a = \k)(l\ + \l)(k\ with a = 1, 15 according 
to the ordering of the pair. Similarly we define 15 antisymmetric combinations 
X@ = i ■ {\k){l\ — with j3 = 16, ...,30. The remaining five operators are 

diagonal in k and I. The multipolar interaction of ea. (|5-7[) is written as 

H ** = " E E Vfx?xf = - £ Xf • % ■ X, , (5-8) 

(ij) a,0 (ij) 

where X = [X , X 2 , .., A 35 ], and i,j are site indices. The susceptibility matrix in 
the momentum space can be expressed as 

X(q) = (l-Xo-^(q))^-Xo, (5-9) 

where V a @(q) is the Fourier transform of Vf* , and xo is the single-site susceptibility 
matrix. The elements are given by 

(Xo)^ = E PiE * ] : P E {El \ k\X«\l)(l\Xe\k), (5-10) 
k,i 1 k 

where p(Ek) = exp(— /3£'/ c )/^ fc exp(— (3E^). In the right-hand side of ea. (|5TU|) . we 
only need terms with k 7^ I, since the interaction Hamiltonian 1)5 -7|) does not have the 
diagonal part of operators X a . The calculation of (xo) al3 can easily be performed 
with use of the magnetic eigenstates, which are nondegenerate except at the level 
crossing point. 
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In the presence of ferromagnetic interaction, the effective field acting on CEF 
states can be stronger than the external magnetic field. It is then important in 
deriving the phase diagram to include renormalization of the external magnetic field 
h = g^sH. We consider the effective magnetic field given by h e g = h — £<?dip(«^iil) 
where the moment (Jm) along (111) depends on h e g and temperature T. The 
effective field should be determined self-consistently. Actually we first take h e g as a 
given field and derive the matrix xo under h e g. To each value of h e s, we determine 
h by the relation h = h e g + zgdip^i in- 
susceptibility matrix. This is equivalent with the condition 



The solution of eq. (|5-lljl gives the instability condition of the paramagnetic phase. 
Note that the phase boundary may not coincide with the instability line if there is a 
first-order transition. In this paper, we do not go into detailed inspection about the 
order of the transitions. 

5.3. Ordering vector and instability lines in magnetic field 

Under the reduced symmetry with finite magnetic field along the (111) direc- 
tion, the original system allows only two different symmetries for the order pa- 
rameter: the Cs(ri) singlet and the C^r^) doublet representations. 23 ^ Within the 
two-dimensional local Hilbert space, which is formed by the two states near the level 
crossing point at H = H cross , the order parameters with 63(13) symmetry have 
nonzero inter-level matrix elements. Therefore, if the high-field phase has the same 
ordering vector as the low-field AFQ phase, both phases can be connected smoothly. 
If, on the other hand, the high-field phase has a different ordering vector q, the two 
phases cannot be connected smoothly. 

We derive the instability from high temperatures toward the low-field and the 
high-field phases with several sets of interaction parameters. We study both cases 
with the ordering vector q = (1,0,0), which leads to a two-sublattice order on the 
bcc lattice, and q = 0. The quadrupolar coupling constant <? qua d is determined to be 
8.93 mK by the zero field transition temperature Tq = 6.5K. The small value of g^ad 
as compared with Tq is firstly due to the large numerical value of the quadrupole 
moment, and secondly due to the large value of z = 8. Figure |H1 shows examples of 
the calculated phase diagram with the following parameters: W = — 0.9K, x = 0.98 
and y = —0.197, leading to Aqef = 0. Two different types of solution can be 
obtained depending on the magnitude of ferro-type coupling constants: (a) both 
low- and high-field phases have the same ordering vector q = (1,0,0); (b) the high- 
field phase has the ordering vector q = 0. In the case (a) we obtain either continuous 
phase boundary or two separated phases. For example, with g^ p = — 125mK and 
5oct = — 25.7mK we obtain the continuous behavior as in the left panel in FiglHl 
while with the same g^ip but with g oct = — 25.9mK, two phases are separated. In 
the right panel, on the other hand, the high-field phase has the ordering vector 
q = 0. In this case, the the high- and low-field phases are separated by the first- 
order transition. The effects of the ferro-type interactions on the phase boundary 
depend on the multipolar matrix elements and the mixing among multipoles due to 




(5-11) 
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Fig. 8. Temperature-magnetic field phase diagram with two alternatives, q = (111) or q = 0, of the 
ordering vector in the high-field phase. Dashed line shows the observed phase boundary. The 
parameters used are g qua ,d = 8.93mK for both left and right panels, and Left: <?di P = — 125mK 
and g oct = — 25.7mK; Right: (?di P = — 102. 5mK and g oct = — 34.2mK. 

the symmetry lowering as we increase the magnetic field. We found that the dipolar 
matrix elements are effective in reduction of the transition temperature in the low- 
field regime (H < H CIOSS ), while the T 5u octupolar matrix elements have larger 
effects in the high-field regime (H ~ H cross ). This latter behavior arises because the 
inter-level octupolar matrix elements between 1^-/23 levels become significant as 
we increase the magnetic field. In the case of larger ferro-octupolar coupling, the 
high-field phase thus takes q = as the ordering vector. 

5.4. Quadrupole patterns in the ordered phase 

In the presence of AFQ order, the relevant triplet states are split. By applica- 
tion of magnetic field, staggered magnetic moments are induced in addition to the 
homogeneous moment. To analyze the splitting of the triplet, we make use of the 
^3-type quadrupolar moments Of, and 0\ projected to the cubic I4 or I5 states. 
In contrast with the /5-type AFQ realized in PrOs4Sbi2, the /3-type AFQ does 
not involve the singlet T\ as seen from eq. ([3-6j) . Hence we concentrate here on the 
triplet only. We write \r t ,m) with m = ±,0 simply as |m), and represent them 
as {|«i), 1 0:2) , I CK3) } = {|0), |+), |— )}. Similarly the cubic i~4 triplet states are rep- 
resented by \ai) (i = 1,2,3). Then we define the 3x3 matrices O2 and C 1 ^ with 
li = 0, 2 by 

= (ai\0%\a k ), (C$) ik = (Oj|0^|o fc ). (5-12) 
By comparing matrix elements in both basis sets, we obtain 

0° 2 = a Q C Q 2 + b Q C 2 2 , 

0\ = a Q Cl - b Q C° 2 , (5-13) 
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where = (9w— 2)/7, 6q = — y/3w(l — w)/7. The limits w = 1 and give the cubic 
triplets i~4 and I5, respectively. Indeed for w = 1, we recover 2 = C 2 and C| = C| 
from (|5-13[) because of ciq = 1 and 6q = 0. The states |0), (|+) + \—))/y/2 and 
(|+) — |— ))/v / 2 are simultaneous eigenstates of Of and 2 . The eigenvalues of 2 can 
be obtained from the eigenvalues — 14/v3(l, 1, —2) of C 2 , together with eigenvalues 
±14, of C 2 by proper combination. Generally, the eigenvalues of 2 consist of three 
singlets. However, as shown below, the eigenvalues under the quadrupole molecular 
field consist of a doublet and a singlet. 
Using eq. Q5T3JI we obtain 

(Olf + {Olf = (4 + b 2 Q )[(C° 2 ) 2 + (C 2 ) 2 ] = 7 -f{a 2 Q + (5-14) 

which is a scalar matrix in consistency with the fact that (0 2 ) 2 + {O 2 ) 2 is a Casimir 
operator. The factor Oq + bg = (4 — 15d 2 + 60d 4 )/49 reduces to unity with d = 1, 
and to 4/49 with d = 0. Namely the i~5 triplet has a smaller quadrupole moment as 
compared with the I4 triplet. 

In considering the AFQ ordering pattern, it is convenient to begin with a two- 
site problem with sites A and B at zero temperature. The interaction Hamiltonian 
is given within the subspace by 

H(r 3 ) = Ar 3 (0 2 %0 2 % + 0\ A 0\ B ). (5-15) 

We take the nine-fold basis |A;)a|0b f° r the two sites, where k,l can be 0, + or — . 
By diagonalization of 7i{r^) we find that the nine- fold degeneracy splits into a six- 
and a three-fold multiplets. With antiferro-type interaction (Ar 3 > 0), the ground 
state is the six-fold multiplet with states 

|0)^|+)b, |0)a|-)b, \+)a\0)b, Ha\0)b, 
^={\+)a\+)b-\-)a\-)b). (5-16) 

To obtain insight into the nature of degenerate wave functions, we associate the 
triplet states with I = 2 spherical harmonics as 

|0>~Y 2 °~3z 2 -r 2 , |+) ~y 2 2 ~( x + iyf, I-) ~Y 2 2 ~(x-iy) 2 . (5-17) 

If we choose the state |0)a, namely (3z 2 — r 2 )^, the six-fold degeneracy allows com- 
binations \+)b i \—}b to make (3z 2 — r 2 )^(x 2 — u 2 )b and (3z 2 — r 2 )A(xy)B as the 
degenerate ground states. Similarly, combinations of states in the second line of 
ea. (|5-16|) give other degenerate ground states (xy)^(x 2 — v 2 )b and (x 2 — y 2 )A(xy)B- 
We remark that the six-fold degeneracy comes from the number of ways to choose 
different orbitals at two sites. 

With |0), (|+) + \—))/y/2 and (|+) — \—))/y/2 chosen as the basis at each site, 
the mean-field theory at zero temperature becomes exact for the two-site system. 
The absence of off-diagonal elements in the two-site Hamiltonian persists in the bcc 
lattice, which can be separated into A and B sublattices. Hence with only the I3 
intersite interaction, the degenerate ground state of the bcc lattice should be given 
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exactly by the mean-field theory. Recognizing this situation, we consider the mean 
field of W(P 3 ) at site B: 

H B = Ar 3 {(O 2 )aO° 2 + {Ol) A Ol) = Ar 3 {Q A O 2 + g A 2 2 ) , (5-18) 

where we have used the notations (0 2 }a = Qa and (0 2 }a = QA- If we choose the 
state 10)^, Hb is diagonal and has the following structure: 

1 / 2 \ 

h m = r-{(a k \H B \ai)} = 2c -1 , (5-19) 
Ar 3 V - 1 / 

where c = (l4/\/3) (ciQ+&g) and we have used Qa = 28/v / 3a<2 and qa = —28/\^3bq. 
We note in eq. (j5-19j) that the splitting into a singlet and a doublet is due to the mix- 
ing of moments 2 and 2 in the quadrupolar mean field. This quadrupolar field 
can be expressed as 

QaO° 2 + q A 0\ = 28/VH(a 2 Q + b 2 Q )C° 2 ~ 3z 2 - r 2 . (5-20) 

The states |±)b constitute the degenerate ground state. If we choose the states 
(|+) ± |— ))a, on the other hand, the resultant quadrupolar fields are now given by 

Q a O% + q A 0\ = -U/V3(a 2 Q + b 2 Q ){C° 2 ± V3C|) ~ 3y 2 - r 2 , 3x 2 - r 2 , (5-21) 

Thus the degenerate ground state is formed by the states 1 0) b and one of ( | \ — ) )b ■ 
From the consideration above, we conclude that a model with nearest-neighbor 
AFQ interaction of i~3 quadrupoles has a degeneracy with respect to different kinds 
of quadrupolar patterns. In this sense, the I3 antiferro-quadrupolar ordering model 
within the triplet state is similar to the three-state antiferromagnetic (AFM) Potts 
model, 34 ) which possesses a macroscopic degeneracy in the ground state. The degen- 
eracy in our case should be broken by dipolar and/or octupolar intersite interactions. 

Experimentally, the distortion of ligands observed by X-ray or polarized neutron 
diffraction is [5, 8, 5 ]-type. 35 )' 36 ) The dominance of this kind of quadrupolar moment 
is also indicated by the results of polarized neutron scattering 37 ) where the induced 
AFM moments is about 2 times larger for the field direction (0, 0, 1) than for (1, 1, 0). 
Taking 3z 2 — r 2 -type quadrupolar field on one sublattice, however, the only way to 
keep the the local Th symmetry at the other sublattice is to break the time-reversal 
symmetry. Namely, if either |+) or |— ) is realized, it is possible that the lattice 
distortion with an 2 component is absent. However, a large magnetic moment 
should emerge in this case, in contradiction to the experimental results. If the time- 
reversal symmetry is preserved, the only possible combinations are (|+) + |— )) and 
(|+) — |— ))• In this case, a non-zero x 2 — y 2 -type quadrupolar field should emerge on 
one of the sublattices. Namely, the AFQ order should always be a mixture of 2 and 
0\ components in the triplet subspace. It is an interesting future problem to deter- 
mine the accurate AFQ pattern in PrFe4Pi2 by detailed comparison of experiment 
and theory. 
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§6. Summary 

We have studied two elements which seem essential in understanding physics of 
Pr skutterudites: on-site hybridization and intersite interaction. The hybridization 
between 4/ electrons and the surrounding ligands plays a decisive role for the char- 
acteristic CEF splitting and the Kondo-like behavior. Emergence of the Kondo effect 
in skutterudites depends on the nature of the triplet wave functions; whether they 
are dominantly composed of i~4. With this observation we have explained why the 
skutterudite PrFe4Pi2 shows the Kondo effect, while another skutterudite PrOs4Sbi2 
does not. More detailed account has been given in refs.13), 14). 

As concerns typical phenomena brought about by intersite interactions, we have 
proposed a comprehensive picture for the i~3-type quadrupole order in PrFe4Pi2- 
The low- and the high-field phases have been reproduced simultaneously with the 
assumption of vanishingly small splitting between the singlet and triplet states. The 
key point for appearance of the high-field phase is that a level crossing occurs only for 
the magnetic field direction (111) in the pseudo-quartet CEF scheme. The magnitude 
of the dipole and octupolar interactions controls whether the high-field phase has 
the same (q = (1,0,0)) or different (q = 0) ordering vector as that in the low-field 
phase. In the q = case, we expect macroscopic lattice distortion in the high-field 
phase, and the corresponding elastic anomaly. It is highly desired that experiment 
in the near future will determine the type of ordering in the high-field phase. 
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